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SUMMARY 


Calculations are presented for the forces on a thin supersonic wing 
underneath which the air is heated. The analysis is limited principally 
to linearized theory but nonlinear effects are considered. It is shown 
that significant advantages to external heating would exist if the heat 
were added well below and ahead of the wing. 


INTRODUCTION 


Efficient propulsion by a ram-jet engine which brings the air to a 
low subsonic speed for ignition and burning is, at present, limited to 
the Mach numbers below which the stagnation temperatures do not present 
prohibitive cooling problems. One method of propulsion that diminishes 
the cooling problem and still makes use of the surrounding air is to add 
heat directly to the external supersonic flow. The object of this paper 
is to study the latter effect and to find heat distributions which are 
in some sense optimum. Specifically, calculations are presented for the 
forces on a thin supersonic wing underneath which the air is heated. The 
analysis is based on linear theory so the flow disturbances must be small. 
However, the importance of certain nonlinear effects are discussed. 

Other authors (see refs. 1, 2, and 3) have analyzed exterior heating 
by means of linear theory. The conclusion reached in these papers, as 
well as in reference 4, is that no substantial gain in airplane range 
performance can be expected from small amounts of heating in a thin (on 
the order of the wing thickness) layer beneath the wing at supersonic 
(H» = 3 to 5) speeds. The analysis presented herein is in basic agreement 
with such a conclusion but it shows that this kind of heating, while per- 
haps the most practical at present, is far from optimum. It is true that 
the optimum heating in a supersonic stream may be attainable only by some 
novel form of energy release such as electromagnetic radiation, but the 
attractiveness of external heat addition appears to be limited only by 
the extent to which it can be controlled. 
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speed of sound 
wing chord 

specific heat at constant pressure 

specific heat at constant volume 
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drag coefficient, 
lift coefficient, 
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pressure coefficient , 
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net force coefficient in free-stream direction 
drag 

integrals of heating function (see eqs. (52) and (56)) 
maximum depth of heated region 
see equation (76) 
lift 

Mach number 
local pressure 
power 

power coefficient. 
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dynamic pressure, — pU 2 

heat added per unit volume per init time 
dimensionless heat term (see eq. (27)) 
length of heated region measurec. from trailing edge 
gas constant 


wing thickness 
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a,Q 


u,Z 


c,t 


stagnation temperature 

thrust from added power 

total velocity in the x direction 

perturbation velocities in x,y,z directions 

dimensionless wing volume, — 1 ™ ie ^ (unit span) 

c 2 

total velocity (vector) 

Cartesian coordinates 

distance from z = 0 plane to wing surface 
wing angle of attack 

Ao 2 - 1 


ratio of specific heats 
engine efficiency factor, 




engine 


wing slope 

heat added per unit air mass 
local density 
wing thickness ratio 
perturbation velocity potential 


Subscripts 

free-stream conditions 

conditions at initial and final heating stage 
velocities caused by aerodynamic and heat effects 
conditions on upper and lower surface of wing 
slopes attributed to wing camber and twist 
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BASIC ASSUMPTIONS AND EQUATIONS 
Equations for Small Perturbation Flow 


Neglecting viscosity and heat conduction, and assuming a perfect 
gas, we can write the equations of state and continuity of mass, momentum, 
and energy in the form 


p = pRT 

(la) 

(V • grad)p + p div V = 0 

(lb) 

p(V • grad)\f + grad p = 0 

(lc) 

. grad) (c p T + ^ 2 = p 

(Id) 

added per unit volume per unit time. 

If the flow 


is limited to one having only small disturbances, these can be combined 
to form (see refs. 5 and- 6) 

= - Q V V7 " 1)M “ 2 




( 2 ) 


where cp is the perturbation velocity poteitial and is the free- 

stream dynamic pressure. The term Q v is related to p, the heat added 
per unit mass of air, by the equation 

(V • grad)p = (3a) 

which, for small disturbances, yields 

n = 7-Tr- / QvUi^zJdXj. (3b) 

p JWg(y, Z ) 

where x = g(y, z) is the equation of the forward surface of the heated 
region. Similarly, the equations for press are and density reduce to 


P-Poc _ 2u 
%> " U 00 
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rp + (r-i) 
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( 5 ) 
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The basic assumptions are that the effects of viscosity and heat 
conduction are negligible and that the air is a perfect gas. The mechanics 
of the heating itself are ignored and assumed to be independent of the gas 
behavior. By far the most restrictive assumptions are those which limit 
the flow to small disturbances. Since we are interested here only in high 
supersonic flows , these can be expressed as 

« 1 ( 6 ) 


(y-i)d 

aoo 2 


« 


1 


(7) 


where t is a representative slope of an exposed surface. The first of 
these conditions is the usual one that limits the application of linear 
theory when applied to high Mach number flows and the second one represents 
the limitation on the heat intensity in such flows . 


Discussion of Validity of Linearized Equations 


The conditions (6) and (7) which have been imposed on flow fields to 
which equations (2) through (5) apply, have a clear interpretation; namely, 
that the wing slopes be small compared to l/^ and the heat addition be 
limited to 


7-i m^ 2 r x 

Q v dXl<<1 (8) 

§ 

Quite possibly other special types of heating not restricted by (8) 
are governed by the linearized equations. Consider, for example, a one- 
dimensional flow in which there is no area change. In this case the exact 
equation (neglecting viscosity 


aoo 2 + (y-l) ( q + 


which has the solution 


and heat conduction) can be written 


U 2 -U 2 ' 


- ir 


fU = ( 7 _i)u p 

dx dx 


(9) 
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( 10 ) 


if u-0 at U = Uoq. In terms of the perturbation velocity u, the 
pressure and heat input can be expressed as 


P-Poo _ 2u 
3oo " 


( 11 ) 
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and 


7- 1 

2 





x 2 

7+1 f _u_) 

Uoo/ 


From the linearized equations, 


xx 


(Hx» 2 -D<P ; 

P-Pqo 

u m 


on the other hand, 

(7-DQ^ 2 

2u 


one finds 


( 12 ) 


(13) 

(14) 


and 


7-1 

2 




(15) 


Clearly, in the case of one -dimensional flow, linearized theory 
predicts the heating effect on the pressure to a higher order than that 
given by expressions (6) and (7). In the oie -dimensional case the linear 
theory is valid for supersonic Mach numbers if u/Uoq « 1 and not 
H»u/Uoo « 1. We will return to this point presently. 

An obvious way to check a proposed approximate theory is to apply it 
to examples which have been studied by more exact methods. One such set 
of examples is supplied in reference 4 wherein the effect of adding heat 
under a two-dimensional, 5-percent -thick, biconvex wing at 0° and 2° angles 
of attack is calculated by a nonlinear theory. The heat was added so that 
the stagnation temperature rose linearly from its free stream to a speci- 
fied final value. Figure 1 gives a comparison of these results with those 
found by the linear theory 1 used in this report. 

The two theories are in good agreement especially as to the effect 
of heat addition on the pressure. This agreement is surprising in the 
case of figure l(a). One can show that 


(/-l)h 

o 2 




:>x i 


X -P -X ! 



(16) 


and for x = x f , = 5, and Tf/Ti = 1.243 this becomes 1.458 which 
certainly violates condition (7). Yet in figure l(a) the linear theory 
follows the more exact one even to the trailing edge. 


1 The analytical form of the linearized results and a brief discussion 

of their theoretical development are given in the next section. 
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Stream tines in isen tropic 


Perhaps one explanation for this degree of accuracy can be drawn 
from the discussion of one -dimensional flow. It is well known that, to 
a first order, the stream-tube area intercepted by a Mach line does not 
change in the region bounded by the fore and aft shocks in a two- 
dimensional flow, as shown in sketch (a). The linearized treatment of 
the heated flow contains the result 
that under the condition of constant 
heat addition, as is the case for both 
examples in figure 1 , the stream tubes 
are unaffected by the heating in 
region A of sketch (a), and begin to 
diverge at a rate proportional to the 
distance from the wing in regions B 
and C. The heat added in region C, 
however, has no influence on the wing 
pressure distribution and the value of 
( 7 - l)n/aoo 2 at p, the last point on the 
outer portion of the heated region 
that can affect the wing pressures, is 
about 0.66. Thus, over most of the 
part of the heated region that contri- 
butes to the wing pressure field 
either the flow is nearly one- 
dimensional or the magnitude of 
(/-lJli/a^o 2 is small. 



Region of heat addition 
in linearized theory 


* 

B 


From another point of view, the Sketch (a) 

high degree of accuracy for the line- 
arized results shown in figure 1(a) must be considered fortuitous. Mager 
in reference 1 , for example, uses the identical linearized theory and cal- 
culates pressures about 40 percent higher than those shown here. The 
reason can be explained with the aid of sketch (b). The upper part of 
the sketch shows the exact boundaries of the 
wing and heated region. The original height 
h Q of the stream tube which enters the heated 


region with height hi is, in this case, about 
40 percent greater than hi . Mager cast his 
equation for pressure in terms of the power 
coefficient given in reference 4 which, in turn, 
is correctly based on h Q . Hence his calcula- 
tions show the correct value for the power but 
high values of pressure. In this report the 
equations were derived in terms of the boundary 
conditions in the lower part of sketch (b) and 
predict the correct value for the pressure but 
a power 40 percent low. 

The point of all this discussion is that 
figure 1 (a) suggests an accuracy for the line- 
arized theory of heat addition at Mach n umb ers 
and heat intensities which undoubtedly exists 
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only for special cases. In the following sections some quite general 
results concerning the forces developed on two-dimensional wings in heated 
flows are derived. In general their validity depends on the extent to 
which expressions (6) and (7) are true. 


Solutions to the Linear Equations in Two 
and Three Dimensions 


Obviously equation (2) is either elliptic or hyperbolic depending on 
whether is less or greater than 1. By the application of Green's 

theorem, solutions can be written for either the subsonic or supersonic 
case. We will consider here only supersonic flow for which the potential 
can be expressed in the form 


cp(x,y,z) = 


S 


Sip 

Sv 


- cp 


So 

Sv 


dS 


_1_ S 

2« Sx 


JF° 


(7-DV^v 

2c L«> 


dV 


(17) 


where 


Sep 

Sv 


Sep 

Sx 




(18) 


and the direction cosines v ± , v 2 , v 3 of the conormal v are related 
to those of the inwardly directed normal by 


— n^3^ — Av ^ y n 2 — Av 2 * ng — Av 3 


( 19 a) 


A is fixed by the condition 


v i 2+v 2 2+v 3 2 = 

and a is given by the equation 


( 19 b) 


a = arc cosh 


X -<1 


y( y-y x ) Mz-zJ* 


( 20 ) 


The symbol V represents the volume enclosed by the surface S over 
which cp and Scp/Sv are integrated. Applications of equation ( 17 ) to 
supersonic wing theory when Q v is zero are well known. 


In two dimensions the solution simplifies to 


f A t |ds| *11 


(r-i) m/Qv 
2%> 


dS 


0 


( 21 ) 
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where S is now the area enclosed within the curve s. Notice that 
when the integration is along a streamline A = 1, and when along a 
characteristic A = 0. 

In both two and three dimensions the effect of heat addition on 
small disturbance flow fields is identical to that of distributing simple 
fluid sources of strength (7-l)H X) 2 Q v /2q O0 throughout the flow. (This 
effect was pointed out in both refs. 5 and 6.) 


APPLICATIONS 

Heating Under a Biconvex Wing 


Consider the wing shown in 
sketch (c). Its upper surface is flat 
and its lower surface is given by the 
equation 

Z = - ^ x(c-x) - cxx (22) 

where a is the angle of attack. If 
there is no heating, the linear! zed- 
theory value of the pressure coeffi- 
cient on the lower surface is 



or 



(23) 

(24) 


If heat is added in the region shown in the sketch in such a way 
that h, the stagnation enthalpy, increases linearly with x then 


p = Zh 


_ x-x n - 
P x f -x ± 


(Tf-Ti) 



(25) 


where and Tp are stagnation temperatures at x^ and x^, respectively. 
From equation (25) one finds 


p-ph/qv f 1 + i± m A (1 t . A 
2^ A 2 ^ y x f- x i 


(26) 


that is to say, Q, v is a constant over the entire heated region. 
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As an example of the way in which 
equation (21) can be used, consider the prob- 
lem of finding the pressure at the point p 
indicated in sketch (d). Introducing the 
dimensionless heat parameter where 


d y 
6 \ 


( 7 -l)t^ 2 h f Q v 


(27) 


we have 


Sketch (d) 

-£ l as + h^ axi d P l a3 + fe JT_ 


Q /V .dJi 1 dz 1 = 0 


r a scp _ r d dq> . 

4 as J. as 


Ph f jj -- 1-1 

paoep 




padep 


since cp is a constant along oe. These equations give 

\X-(3z 


29 p - ^P a " p 


dx x + 

Szj Phf 


2CPm - 2f P a + 


U„ 


a T phf 


JJ" Qlv&x-i&Zi = 0 

aoep 

JJ' Q, v dyL 1 dz 1 = 0 


padep 


Eliminating cp a , we find 
r x-(3z 


Uqo 

2gh f 


(28) 


(29) 


( j(T 5 ^ dz x + Qv^id-Zi 

\jpabfp abca , 


(30) 


The first term in equation (30) gives the effect of the wing on the 
flow and the second term that of the heat. If Q v is a constant, as in 
our example, the effect of the heat is projortional to the sum of the two 
areas pabfp and abca. The pertinence of these areas is at once obvious 
if we think of solving the problem by using an image system of sources 
having the wing as a plane of symmetry. 

The pressure coefficient on the lower wing surface is given by 
equation (4) and, using equation ( 26 ) for Q v , is 


C 


P 




7x-x±\ 
Xf-Xi ) 


(31) 
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This equation is valid for x± < x < x± + phf , that is until the forward 
Mach line from p intersects the lower front comer of the heated region. 
For + (3hf < x < xf , one can show 



and similar results can be written for x > xp. 

In figure 1 the above results from the linearized theory are compared 
with those from the more exact theory given in reference 4. The close 
agreement between the two theories has been discussed in a preceding 
section. 

We can estimate the significance of the result of external heating 
by studying its effect on airplane range performance. From the form of 
the Breguet range equation 


Range 


f l\ f : fuel heat energy\ , weight full 
y P / \ D / \P er unit weight J weight empty 


( 33 ) 


let us consider the combined term 



where 

T^ total thrust 

D total drag 

L total lift 

P total power input 


In the following, 
when (LUoo/P) d=t ^ 


we make the assumption that maximum range is obtained 
is a maximum and will, therefore, refer to this quantity 


as the range efficiency factor. 


As a simple example let the wing, again, as in sketch (c), be flat 
on top and biconvex on the bottom. For simplicity, however, Q v is now 
taken to be a constant over the rear half of the lower surface in a region 
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Sketch (e) 


hounded fore and aft by Mach lines as shown 
in sketch (e). The pressure on the lower 
surface of this wing is given by equation ( 32 ), 
from which the sum of the forces in the free- 
x stream direction 2 on the upper and lower 
surfaces is 


Cx = Cd o + if t2 + l a2 ‘ | Q v(2t-cc) 


and Cp is the net airplane drag exclusive of the wing wave drag. The 
lift is assumed to come entirely from the wing angle of attack and the 
heating, so 


L __ 4c 

<*oo c ~ 3 


4a 1 tt 
= 3 + 3 


Finally the power is given by 


P = // Qydx <iz 


P 


(7-l)H»" 


Let us now evaluate the performance of the wing and external heating 
in terms of the parameter (LUoo/P)^ _q where C x = 0 implies that the total 

thrust and drag are equal. First it is clear that Cj, and P are independ- 


ent of t, so t should be chosen to minimise 
result 

3Q V 

T - 32 


in equation (35). As a 


(3C X = pen + ba 2 + ^ 


2 In keeping with linearized theory the angle between the normal and 
the vertical is neglect ed in evaluating the lift and drag. 
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For a given power, that is, a given value of equation (40) can be 

solved for a under the condition that C x = 0. For optimum (LUoo/P) 
one finds C x=° 



The results for (LUoo/P)^ _q are summarized in sketch (f) under the 


assumption that The maximum value of (LUqo/P) 


obtained 






Sketch (f) 

by this particul ar me thod is seen, from equation (43), to be 
C(7-1)H»72P](1 -h^/ 2) or about 0.5 Moo. This forms the left-hand boundary 



in the sketch and represents either an infinite heating rate or a zero 
value of C D . The right-hand boundary can rot be exceeded without 

providing some extra thrust by an engine to satisfy the requirement 
C x = 0. 


Arbitrary Heating in Two-Dimensional Fields 


The preceding section contains an analysis of a specific wing section 
heated in a specific way. There is no reason to suppose, however, that 
the particular section and type of heating studied were optimum from the 
point of view of, say, fixed volume and maximum (LUqo/pJq^q. We will show, 

in fact, they are not. In order to prove this, let us now derive the 
simple expression for the lift and drag on arbitrary, planar, two- 
dimensional wings in fields with arbitrary heat addition. 


Evaluating planar wing velocity fields by heat images . - It has been 
determined that a region of heat addition disturbs a linearized flow in 

exactly the same manner as if the field were 
filled with fluic sources of strength 
(Uoo/hp )Q V . Thus if, as in the upper half of 
sketch (g), region A is a region of heat 
addition, the value of u and w at P would 
be 




" = 7 ‘kf On \> 41 


U = :: — 


Sketch (g) 


The l/2 results from the fact that the source 
strength is proportional to l/2 the stream 
deflection in the upgoing or downgoing curves. 

Now place a wing in the plane BB above 
the heated region A. The boundary condition, 
insofar as the heat is concerned, is simply 
that the value of w along BB is zero; the 
distortion of the field caused by the wing T s 
own shape can be superimposed on this solu- 
tion for the complete determination of the 
disturbance fielc.. It is clear then, that 
the effect of the heat addition beneath a 
wing is found by placing an image of A 
above the wing, if in sketch (g), and dis- 
carding the wing Obviously, the pressures 
in waves 1, 2, and 3 in the sketch are the 
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same as if a single region of heat addition were present, hut the 
pressure in region 4 is the same as that for the two heated areas. 

Forces on wing obtained from surface pressure distributions .- Using 
the concepts developed in the preceding section, let us introduce the 
following definitions : 

u a , w a the velocities induced by the wing with no heating in the field 

uq, wq the velocities induced beneath the wing by the heated region 

and its image above the wing without the wing being present 


Then assuming that the heat is 
applied entirely within the fore- 
waves from the leading and trailing 
edges as in sketch (h), one can 
write the expressions for lift, 
power, and drag developed beneath 
the wing in the form (on the wing 
u a = -W a /0 



Sketch (h) 

if = - 2 / ( u Q+ u a ) dx 


P_ 

%> 


~2p 

(r-Dn/ 



^r = 2p X (Ua2+u&UQ)dx 


(44) 

(45) 

(46) 


These expressions will be used in a subsequent section to derive 
optimum airfoil shapes for given heat fields. 

Forces obtained from momentum flux .- Although equations (44) through 
(46) are sufficient to analyze the performance of a two-dimensional wing 
in the presence of heat, it is instructive to re-examine the derivation 
of the drag on the basis of the momentum crossing an enclosing surface. 
This derivation is much more involved mathematically but it serves at 
least the following two purposes: one, it illustrates the existence of a 
new form of "momentum force" in heated flows appearing in addition to the 
usual wave drag and being measured in the wake of the heated region; and 
the other, it is capable of generalization to three dimensions. 



1 6 


Consider the system illustrated in sketch (i). The control surface 
is represented by the dashed lines. Linear theory yields no disturbances 



along the lines counterclockwise from b' to a and from a' to b. The 
momentum drag in the wave from a to a' is well known in its two- 
dimensional linearized form to be simply 

U< - - ^ y ~ = 2p J (u a +u<,) 2 dx (47) 

cl 

A comparison of equations (47) and (46) immediately shows that another 
force having a similar order of magnitude must exist in the system. This 
force is caused hy the change in density (see eq. (5)) behind a heated 
region. 

One can show by retaining the lowest orc.er magnitude terms contrib- 
uting to the drag that the linearized value c»f the wake drag in three 
dimension is given by 


^wake ” P a 


(Moo 2 - 1 )u 2 +V 2 +W 2 


j + py-l)u l i^ ,j | ( Zz±)W» J X u 


d|i 

6x x 


dx n 


r dS 


(48) 


where the dS integration is made over an entire yz plane infinitely 
far behind the wing. In a two-dimensional flow the induced velocities 
behind the trailing-edge wave vanish so the first two bracketed terms in 
equation (48) vanish. The third term, however, depends upon the disturb- 
ance in the heated region itself. Hence, in two dimensions 

b 1 x 

Dwake = Poo^" 1 ) fe f dz i f u |r ^ 

u b -oo i 


or 
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= h juzJJ (u Q +u a )Q v dx i dz i (49) 

s 

where S is the area of heat addition shown in sketch (i) (and does not* 
include the heat image area). 

The total drag is the sum of equations (47) and (49)- This total 
drag must he zero when the wing is a flat plate at zero angle of attack, 
that is, when u a = 0. Hence 


and 



(50) 


The value of u a beneath the 
wing is constant along the 
downgoing characteristics. 

If, in the second term on the 
right-hand side of equa- 
tion (50), we change the coor- 
dinate system so that one of 
the coordinates is parallel to 
a characteristic line and the 
free stream is the direction 
of the other, Q v can be inte- 
grated along the characteris- 
tics. Using equation (30), 
one can show 




S 


u a Q v dXidZi 



1 


and the final expression for the drag is 
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= 2p J u a 2 dx 1 +4p 

where the lengths cc' and dd' are defined in sketch (j). 

Equation (5l) gives the same value of drag as equation (46) since 
Uq along cc' is half its magnitude along the wing as illustrated in 
sketch ( j ) . 


J 


U a UQdX] 


(51) 


Optimum Wings in Heated Supersonic Fields 



Consider the unit chord wing 
and heat field shown in sketch (k). 
The heat is added entirely within 
the fore waves from the wing lead- 
ing and trailing edges since heat 
added outside this region will not 
influence the wing and would not, 
therefore, exist in an optimum com- 
bination. Consider, also, only 
cases where heat is added beneath 
the wing. Then equations (44) 
through (46) can be expressed in 
the form 


C - 2(3 
C D - — 


U, 


oo O 
1 


u z 2 (x) + U u 2 (x) - ^u z (x)G(x) 


&X 


°L " tF 

u oo ■ 


u z (x) - u u (x) - ^ G(x) 


dx 


P 


c 


? - r 

b-DM^J o 


l 

G(x)dx 


where 


G(x) = J Q^(x-p :,z)dz ( 52 ) 

and a and b are the extremities of the heated region along the charac- 
teristic line passing through x on the airfoil. Introducing the fact 
that u z /U OT = A z /P and u^U*, = -A u /p where A z and A u are the lower and 
upper airfoil slopes and introducing the notation 
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\a = ^c+H 
^Z = ^c"^t 


(53) 


A c being the slope of the airfoil camber line and At being the upper 
surface slope of the symmetrical thickness distribution, one can show 


pc x = 


3C Do +2 / [2A c +2A t -(A c -A t )G(x) ]dx 


(54a) 


Ct = - 7T 


P J L ~' C 

O 


[ 2A C -G ( x ) ] dx 


(54b) 


P c = 2 / G(x)dx 


(54c) 


The coefficient 


is the net force in the stream direction and 


is positive when a net drag results. The coefficient Cp o is the net 
drag to be overcome by the heat introduced under the wing exclusive of 
the wing wave drag. In the absence of any other form of propulsion it 
represents the friction drag, induced drag, and the wave drag of the 
other airplane parts. Our object is again to seek a maximum value of 
LUoq/P when C x = 0. 

Since and P are not affected by A-fc* we can start by finding 

the thickness that minimizes C x . If we require a net volume V (dimen- 
sionless, equal to (volume/c 2 times unit span)) and insist that the 
airfoil close, this leads to the shape 

A t (x) = | (6V-3G 1 + 2G 0 ^) - ^4V-2G 1 + G 0 ^ - G(x) (55) 

where 


G n = / G(x)dx 


( 56 a) 


G 1 = J' x G(x)dx 


( 56 b) 


Then if A c is optimized by minimizing the quantity 


I = C x + AC L 
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where A is a constant, there results 





(57) 


and. 


pc x = 3C Do +A 2 +6v( 2V+G 0 -2G 1 )4G 0 2 -3G 0 G 1 +3G 1 2 - G 2 (x)dx ( 58 a) 

°o 

C L = | (G q - 2A) (58b) 

P c = — 2r -~— (58c) 

(7-1)Hx> 


If heating is absent, equation (55) gives the familiar result that 
the closed airfoil having a minimum drag for a fixed volume has a biconvex 
section. The optimum section shape with heating depends, obviously, on 
the nature of the heat addition, but even if G(x) is constant over a 
given portion, the airfoil is, in general, e combination of a double-wedge 
and biconvex section. 

In the interesting case when V it sell is chosen to minimize C x 
(under the supposition, which must be verified, that the resulting section 
is real), there results 


v = - 1 G ° + 1 G 1 

(59) 

r x 


Zt=r G o-r g( j 1 )dx 1 

( 60 ) 

J o 


r x 


z c = 1 A + % J G(X: )dXl 

n 

(6i) 


o 


where z = Z^(x) and z = Z c (x) are the equal ions of the upper surface of 
the thickness distribution and of the camber line, respectively. Obviously 
a real airfoil exists if 


G(x)dx > 



o 


x 


o 


( 62 ) 
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for 0 < x < 1 . The equations for the drag, lift, and power become 


PC X = |3C D + A 2 - | / G 2 (x)dx + i Go 2 

(63a) 

* J o 4 


Cl = | (Go - 2A) 

(63b) 

r 2 Go 

r C - 2 

(63c) 

( 7-1 )H» 



and, further, the range efficiency factor becomes 



External Heating Only 


Let us study the case when is a consta nt over the area 

illustrated in sketch (z). Then setting 7 Hx > 2 -1 ~ Hx> 



An inspection of equation (6 5 ) yields the 
principal message of linear theory with 
regard to heating in a supersonic stream. 
Notice that no matter how small r 
becomes Q y can be made large enough 
(and to. stay within the bounds of linear 
theory, this means that the depth of the 
region can be made long enough) so that 
rQ v is fixed. The term 2 /r is then 
dominant and the theory predicts indefi- 
nitely large values of (LUoo/P)^ 

The physical reason for this is 
quite apparent. Two-dimensional linear 
theory predicts no attenuation of 



Sketch (z) 
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pressure along Mach waves. The efficiency of a unit heat source in 
producing pressure on the bottom surface of the wing is, therefore, inde- 
pendent of the source position so long as its upward going Mach wave 
strikes the wing. In the prediction of an optimum, then, the theory will 
lead to the extreme condition of placing all the heat source on the Mach 
line intersecting the point of maximum slope on the lower surface. Obvi- 
ously this result is absurd from a practical standpoint, since the linear 
model would no longer give a valid estimate of the physical flow from 
several points of view. It gives, however, the trend of the first-order 
terms in forming a criterion for optimum heating condition in supersonic 
flow and provides a good starting point for the study of higher order 
effects. A detailed study of a particular case is presented in a 
subsequent section. 

Let us look at equation (65) in a different light and find the value 
of r which requires the minimum heat input for a given (LU ot /P)q x _q. It 

is apparent from sketch (m) that these minima exist. Considering this r 



Moo 



Sketch (m) 

to be optimum, we can construct the curves shown in sketch (n) where Xj 
and X 7 are the slopes of the front and rear portions of the lower sur- 

face, respectively. Under the assumption that Qy is a constant over 
the region shown in sketch ( l ) , these curves give the lift and power 
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Sketch (n) 


coefficients, the extent of heating, and the surface slopes for the 
condition of minimum heating required to produce a given (LU^/P)^ 
at a given Mach number. x 


Engine Combined With External Heating 


If we seek to find the effect of combining engine thrust with the 
lift and thrust provided by external heating, we can return to equa- 
tions (5*0 in which C x is not zero and the power coefficient is 
increased by C x / r\ where 




?hUo 


P 


engine 



2b 


is the engine alone efficiency parameter. Combining equation (55) with 
( 5 ^a) and optimizing the volume, one can write 


0 C 3 


1 l 

= 0C Do +2 / [ 2A C 2 - A C G ( x ) ] dx + i G 0 2 - ^ / G 2 (x)dx ( 67 a) 


PC L = -2 / [2A c -G(x) ] dx 


P n = ^ — 2 / G(x)dx + i C x 


(67b) 

( 67 c) 


For simplicity, again assume Q v is a constant over the area shown in 
sketch (z), and, further, that 


Ac 


- 0-0 > 


-a 


1 9 


0 < x <: 1 - r 

1 - r < x < 1 


where a Q and a x are constants. Now optimizing 


( 68 ) 


LUqq _ Ua 0 (l-r )+4a x r-2rQ v 

P 2(3 jy [Cx 

( 7 -l )^ 2 11 


for a fixed 


PC X = pCp o +4a c 2 (l-r)+4a 1 2 r-i-2Q v cc : r - i Q v 2 r(l-r) 


gives 



r Q v + N /r(2- r)Q v 2 +4p((! x -C Do ) 


2prQ v 

(7-DH» 2 



(69) 


When C x , which now represents the thn.st which must be provided by 
the engine, is zero, this reduces to equation ( 65 ). For zero external 
heating it becomes 


( LOO. _ 2 yp( c x- c i >o) 

V p ) 

v 'o 


P^x 

1 


(70) 
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which is a maximum for C x = 2Cp o ; being, in that case 



If we seek the value of C x which makes LUqo/P in equation ( 69 ) a 
maximum, there results for p » 



where 



and (LUqo/P^ has been defined in equation ( 71 ) as the maximum range 
efficiency factor for zero external heat addition. 


Sketch (o) illustrates the general 
nature of the results. If 
(5/H»)(LU 0o /P) o >1 (for an aerodynamic 
lift -drag ratio of 6 and an engine effi- 
ciency of 0.3, (LU 00 /P) 0 = 1.8) a small 
amount of external heating decreases 
the range efficiency factor. However, 
the efficiency can be increased if 
enough heat is supplied so that 

> it r jl _! - 

r(2-r) . 



(75) L 

If (5/Hx,)(LUoo/P) 0 < 1, LUoo/P can be 2 \j C D 0 M *> 

increased by even the smallest amount 
of external heat addition. 


Sketch (o) 



2 6 


Let us set 



V p / 

v / o 


and seek the value of r which requires the minimum external heat to 
produce a given K. We find for (5 K/Moo) (LU 00 /P) o > 1 



(78c) 


Typical results arising from the application of equations (77) and ( 78 ) 
are shown in sketch (p). The results apply to the case in which the range 
efficiency parameter is to be increased 40 percent above its maximum value 
obtainable with no external heating. Two pc ssibilities are considered; 
one, no engine, and the other, a combinatior. of engine and external heat. 

In both cases r, the extent of heating, is chosen to minimize Q v . For 
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values of l^ 3o /(LU 00 /P) 0 below_ 
about 3.46 the magnitude of Q v 
is lover if no engine at all is 
used. Above 3.46 less external 
heat is required if an engine is 
used and the percent of total 
power required of the engine is 
shown in the lower part of the 
sketch. 


Application to a Particular 
Example 

In order to evaluate the 
rather general results presented 
in the preceding sections, let us 
consider in some detail a partic- 
ular example. Let us consider 
the case of external heating only 
and fix at the outset the follow- 
ing parameters 

Moo = 6 

Cn = 0.0034 > 

^O 

(LUco/P ) Cx = 0 = 2 ' 4 

(79) 

The value of Cp Q is not exces- 
sively low when one considers the 
reduction of the turbulent skin- 



friction coefficient with Mach 

number. The value of (LUoo/P)q = q 

x 

represents an equivalent conven- 
tional airplane combination of 
engine T^Uqo/P =0.4 and aerody- 
namic L/D = 6. These values are 
what might be expected 3 at the 

present time of ram-jet airplane performance at 


Mao 



Sketch (p) 

] ^=6. Therefore, if 


the values in equations (79) could be attained for a flying vehicle, they 


would make this kind of heat addition warrant some consideration. 


The values assumed in equations (79) when used in conjunction with 
sketches (m) and (n) represent the following conditions 


3 It is assumed that the boundary layer is for the most part turbulent. 
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r = 0.5 

c L = 0.067 
P c = 0.028 
A u = 0 


A 7 = -0.1 

L f 

A, = 0.1 

L r 








( 80 ) 


The lift and power coefficients are reasonable and the lower surface 
slopes are less than 6° with reference to the free stream. The term 
involving the heating is the really critical factor. 

A principal assumption underlying the linearized theory is that 
(/-l/a^ 2 ) (Sja/Sx) is small. For the value of Q v given above, this means 
that hf must be of the order of a wing chord for the above figures to 
be first-order estimates of the heating effect recorded in sketches (m) 
and (n) under the conditions imposed by equations (79)- The length of 
the heated column is not so surprising if we remember that we are 
attempting to drive an airplane with an unswept 5-percent -thick wing 
through the air at a Mach number of 6 with a:i equivalent L/D of around 
6 and we are attempting to do this by adding only a small amount of heat 
along any given streamline. That such efficiencies are even theoretically 
possible is due, as has been pointed out, to the fact that the heat is 
added to a two-dimensional supersonic stream so that pressure disturbances 
along a characteristic line accumulate with Little loss in intensity. 

A very valid question remains as to whether the low heating rates 
we have considered are optimum even if they are practical. They were 
considered first, not only because linear theory demanded it, but also 
because small disturbance fields are, in conventional aerodynamics, gen- 
erally the most efficient. Quite possibly i; is preferable for the heat 
to be intense enough to produce higher order disturbances; this point is 
discussed briefly in the next section. 


Some Remarks on Nonlinear Effects 


Estimations. - Although linear theory cai be depended upon, for most 
practical purposes, to give reliable values for the disturbances near a 
thin wing if these disturbances are very small relative to the free 
stream, there is always the question as to i;s dependability when the 
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disturbances are those produced by a vehicle actually intended to fly. 

In problems as simple as the ones we have been studying, estimates of 
what are probably the principal errors can be made for particular cases. 
For example, let us examine more closely the true nature of the flow, the 
linearized version of which is described by equations (79) and (80). 

The discussion is divided into three parts; the area of heating, the 
pressures caused by the wing shape, and the pressures caused by the heat 
addition. Quite clearly the very fact that the nonlinear aspects of the 
flow field are being considered means that these three points can not 
really be considered independently; at this point, however, it seems 
reasonable to expect that their independent effects deviate more from 
linear theory than their interdependent effects. 4 This matter can really 
only be settled by solving the exact equations for the given boundary 
conditions . 

Consider now the actual region of heat addition for the wing in our 
example. Because we have applied our boundary conditions in the z = 0 
plane (as is usual in linearized theory) the assumed region of heating is 
that shown in the upper part of sketch (q). In a real flow, however, the 



Linear - theory boundary conditions 



Exact boundary conditions 


Sketch (q) 

width of the heated region required to influence the rear half of the 
wing is reduced to that shown in the lower half of sketch (q). In our 
example r a =j 0.4r. Now according to linear theory, the value of the 
pressure coefficient on the wing due to the heating is 


r ~ _2_ ( ZlL 

p ~ w 2 


h f 


4 Except, perhaps, for the possibility of separation in the heated region. 
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and the power coefficient is 


2 (z± cjA 

<3cc U oo c ( 7 -1 )h/ Vaco 2 ox / 

In other words, the pressures, that is, the lift and thrust, are unaffected 
by shortening the region, but the power required to produce them is 
decreased to 0.4 of the value based on the linearized boundary conditions. 

The pressure coefficients on the bottom of the double-wedge section 
described in equation (80) flying at M = 6 in an unheated flow are, by 
linear theory 

Cp = 0.033, front half 

Cp = -0.033, rear half 

On the other hand, simple wave theory (i.e., Prandtl-Meyer expansion for 
both compression and expansion) gives 

Cp = 0.048, front half 

Cp = -0.024, rear half 

Notice that the drag is increased by the amount 


ACq 


(S|Cpl) 3lrt -(£| c pl>i im: ar 


t = 0.0003 


(and this can be interpreted as reducing the Cp 0 in equation (79) to 
0.0031) and the lift is also increased by the amount 



( ZC p) S wt-< ZC p) 


'P 'linear 


0.012 


Finally, consider the pressure produced by the heat addition. The 
error incurred by the use of linearized theory in this region is difficult 
to evaluate since an exact solution for the particular example under con- 
sideration is not available. A very rough estimate of the error is made 
in the following discussion. 

First, if p is assumed to be independent 5 of V 2 , the pressure 
coefficient can be written 


5 For this example, the effect of p is, in any event, not large. 
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or, on the wing, according to equations (79) and (80) 



where is the value of the pressure coefficient given by- 

linearized theory and x is the distance into the heated region. 

Sketch (r) shows the results. The two trends are clearly indicated: one. 




Sketch (r) 
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at some point on the wing the actual pressure is greater than that given 
by linear theory, and as the linearized value increases (i.e., as the u 
perturbation on the wing increases) this discrepancy increases; the other, 
as the heated region is penetrated beyond such a point (or as its depth 
is decreased) the actual pressure falls relative to its linear value. 

In order to make use of sketch (r), we must estimate the value of u 
along the surface of the wing. If we use linear theory for the heat 
sources but take into account the fact that the "image plane" is tilted, 
there results 


u heat = 

Uoo P 


So, using equation (80 ) one finds 


^linear 


0.1 


heat 


It is impossible to tell what the nonlinear sffects on this calculation 
would be without a much more sophisticated analysis. It is almost certain 
that the value 0.1 is much too low because the entire effect of the heat 
and reflection is to compress the air, and the linear theory is already 
30 percent too low in an unheated flow when it predicts a value of 0.033 
at Nb = 6. It seems conservative, then to estimate (Cp _ ) at 

0.13. 11 


mear heat 


One can easily show that the value of 
by the wing alone yields the expression 


^linear 


a induced on the rear surface 


« - 0,03 


wing 


Hence, by these approximations, Cp 

about 0.1, which, according to sketch (r) 
of about 0.1^ there. 


on the lower rear surface is 
, gives an average value for 


C 


P 


The following table summarizes the above crude corrections to linear 
theory for the airfoil section and heating rates (hf ~ 0.5c) given in 
equations (79) and (80): 
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Pc 

c L 

C Do 

LUqo 

P 

Linear theory 

0.028 

0.067 

0.0034 

2.4 

Plus correction for size 
of heated region 

.0112 

.067 

.0034 

6.0 

Plus correction for actual 
wing pressure 

.0112 

.079 

.0031 

7.0 

Plus correction for actual 
heating pressure 

.0112 

.094 

.0047 

8.4 


Discussion , - According to the table the linear theory is extremely 
conservative for the particular example studied. When the aforementioned 
corrections are applied, not only is the allowable excess drag (Cp 0 ) 
increased by 38 percent but the equivalent L/D (relative to a ram- jet 
efficiency T^Uoo/P = 0.4) is increased to 21. The reliability of these 
calculations has been repeatedly qualified. They clearly indicate, how- 
ever, that this kind of heating, if it could be practically obtained, is 
worthy of study, and really quantitative theoretical estimates must 
include nonlinear effects. 


Cp 

s 


Exact 


As was mentioned at the end of the last section, the optimum heating 
intensity is likely to produce disturbances well above those dealt with 
by linear theory. The reason for this can be demonstrated as follows. 

The effect of the heating is to increase the pressure of the air. In one 
respect adding a given amount of heat to g 
the field corresponds as far as the pres- 
sure field is concerned to the insertion 
of a wedge-shaped airfoil in the flow. 

It is well known that the pressure rise 
per degree for a wedge is greater, the 
greater the angle of stream deflection ^ 

(see sketch (s)). In the usual case of 
aerodynamic design this means, in a gen- 
eral way, that surface slopes should be 
reduced as much as possible to minimize 
the drag. In the case of heat addition, 
however, the opposite is true. Of £ 

course, no net force acts on a heat 
source whose only use is to provide as 
high a pressure as possible. Therefore, 
inasmuch as a heat source is like a 
fluid source the heat rate should not be 
of a low intensity. 



Linear- 


M*6 

00 


8 


Unfortunately the effect shown in 
sketch (s) is not the only nonlinear 
effect caused by increasing the heat 


8 , degrees 

Sketch (s) 
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rate. Another important effect is to increase the local speed of sound 
which, in turn, decreases the pressure. An estimate of this result has 
been presented in sketch (r) and is illustrated by the reduction in Cp 
for a fixed c p linear when x / h f 6 oes from C to 1. Some optimum relation 

between wing-surface and heating-region geometry and heating intensity is 
indicated and it appears probable that the intensity required will be of 
a magnitude to produce higher order disturbances in the flow. 
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(b) Moo = 3> *i/c = 0.356, Xf/c = 0.690 

Figure 1.- Pressures on lower surface of 5 -percent -thick biconvex wing 
with and without heat addition; Tf/Tj_ = o, - 2° y h^/c = 0.07. 
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